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Top Quark Decay at Next-to-Next-to Leading Order in QCD 
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We present the complete calculation of top quark decay width at next-to-next-to-leading order in 
QCD, including next-to-leading electroweak corrections as well as finite bottom quark mass and W 
boson width effects. Especially we also show the first results of the fully differential decay rates for 
top-quark semileptonic decay t — > W + (l + v)b at next-to-next-to-leading order in QCD. Our method 
is based on the understanding of the invariant mass distribution of final-state jet in the singular 
limit from effective field theory. Our result can be used to study arbitrary infrared-safe observables 
of top quark decay with the highest perturbative accuracy. 



Introduction. The top quark is the heaviest fermion 
in the Standard Model (SM), and frequently plays an 
important role in many extensions of the SM. Therefore 
detailed studies of its production and decay are highly 
desirable. Their precise measurements at the LHC will 
be crucial for the understanding of electroweak symme- 
try breaking and also searching for new physics. Due to 
its large mass, the lifetime of top quark is much smaller 
than the typical time scale of hadronization. For this 
reason, the top quark can be treated as a free particle 
in good approximation, and perturbative calculation of 
higher order quantum corrections to its decay rate can 
be performed. 

Within the SM. the ncxt-to-leading order (NLO) QCD 
corrections to top quark decay width, F t , are calculated 
more than 20 years ago [H]. Employing the method 
developed in Ref. Q, the next-to-next-to- leading or- 
der (NNLO) QCD corrections to T t was calculated in 
Ref. [3J, in the limit of m t 3> mw- Later, finite W boson 
mass effect in the NNLO computation is taken into ac- 
count in Refs. [HQ based on the calculations of top quark 
self-energy as expansion in m^r/mf . All the previous cal- 
culations at NNLO concentrate only on the inclusive de- 
cay width, but differential decay rate is also of substantial 
interest, especially when considering the measurement of 
top quark mass p and electroweak (EW) couplings 0- 
In particular, it's an important ingredient in a fully dif- 
ferential calculation of top quark pair production [8j and 
decay at NNLO in QCD. To the best of our knowledge, 
such a calculation has not been finished so far, and is the 
subject of this Letter besides the total decay width of the 
top quark. 

The formalism. We consider the SM top quark decay, 



t 



W + + b + X, 



(1) 



where X represents any other parton in the final state. 
NNLO QCD corrections to this process consists of three 
parts: two- loop virtual contribution (X contains noth- 
ing), one- loop real- virtual contribution (X contains 1 
parton), and tree- level double real contribution (A" con- 
tains 2 partons). While the amplitudes for each part are 



well defined, integrals over the phase space induce in- 
frared singularities, which must be extracted to cancel 
against those from virtual corrections in order to obtain 
a finite result. In particular, the double real contribu- 
tion is the primary obstacle for obtaining fully differen- 
tial NNLO corrections. In the past decade significant ef- 
forts have been devoted to solving this problem, and fully 
differential corrections have been obtained for a number 
of important processes using quite different methods 0- 
IToj . In this Letter, we solve this problem for processes 
of heavy to light decay at NNLO in QCD, using a phase 
space slicing method inspired by a factorization formula 
for heavy to light current in Soft-Collinear Effective The- 
ory (SCET) [ll| . Below we describe our method. 

Cluster all the partons in the final state into a single 
jet. Let t — (pi,+px) 2 /ntf, which measures the invariant 
mass of the jet, the top quark decay width F t can be 
divided into two parts: 



r t = / dr — 



dr-l=T A +T B , (2) 
dr 



which will be treated separately as explained below, tq is 
a dimcnsionless cutoff for t, and r max = (1 — my/ /mt) 2 . 
We set bottom quark mass mt = in the NLO and 
NNLO QCD calculations. Effects of finite mj, are small 
and will be considered later as a correction to the lead- 
ing order (LO) results. In the limit of r — s- 0, only soft 
radiations and (or) radiations collinear to the b quark 
are allowed. In this region, obeys a factorization for- 



mula 
1 



I* 



dr. 



(0) dr 



r 



= H 



xS 



/i dkdm 2 J(m 2 , ii)S(k, n) 



2Ejk 



+ ■ 



(3) 



(0) 



where we have neglected non-singular terms in r. F 
is the top quark decay width at LO, /i is the renormal- 
ization scale, and Ej = (m 2 — m 2 v )/(2mt) is the energy 
of the jet near threshold. H(x,fi) is the hard function, 
which results from integrating out hard modes of QCD 



2 



in matching to SCET. It has been calculated to NNLO 
in a s pjj]. J(m 2 , fj,) is the quark jet function with mass 
to, whose NNLO expression can be found in Ref. [lij ]. 
It can be thought of as the probability of finding a jet 
with invariant mass m, generated by collinear radiations. 
S(k,fj.) is the soft function, which describes the probabil- 
ity of measuring the light-cone component of the momen- 
tum of soft radiations k s -n, where n is a unit light-cone 
vector along the direction of the jet, to be k. It has also 
been calculated to NNLO in Ref. fia ]. 

Using the NNLO results for the hard function, jet 
function, and soft function, we can then calculate Fa 
at NNLO, utilizing Eq. (J3|), up to terms proportional to 
r . For sufficiently small r , they can be safely neglected. 
The most difficult part of double real contributions are 
included in the calculation of the jet function and soft 
function. Note that F^ is infrared finite, because the 
infrared divergences in the jet and soft function cancel 
against those from the hard function. The spin informa- 
tion of the b quark is lost since spin summation has been 
performed in the jet function. But polarization informa- 
tion of the top quark is retained, due to the fact that soft 
radiations do not change spin. In practice, instead of a 
convolution form, it's more convenient to write Eq. ([3]) 
in a product form: 



1 ^ 14, N 

T t dr 

/ TYl 2 

x lim j I d„ + In — i- , u 



(4) 



a, + In ■ 



~1eV 



2E Jf x^J t r(r?) 



where j and s are the Laplace trasformed jet and soft 
function, respectively: 



j In 



dm exp 



eT E TO( 



J(m 2 ,n), 



s In 



dk exp ( — - ? ) S(k,[i), (5) 



and t v /t should be expanded in terms of plus distribu- 
tion: 



L = V)+E ! r 

r rt i n\ 

' n—0 



(0) 



Substituting the NNLO expansion for hard function, jet 
function and soft function into Eq. (j4]) gives a closed form 
solution of dLt/dr at small t. 

Tb is also infrared finite. In fact, 0(a 2 ) contribution 
to it can be obtained from the NLO QCD corrections to 
t — > W + b plus 1 jet, as long as t > 0. In our calculation, 
the one-loop helicity amplitudes for this specific process 
is extracted from the NLO QCD corrections to single 
top production associated with W boson [l6| . The tree- 



level helicity amplitudes are calculated with HELAS [17 1 . 
Infrared divergences in the phase space integral of tree- 
level matrix elements are canceled by adding appropriate 



dipole subtraction terms [l8|. For later convinence, we 
further divided the 0(a 2 s ) contributions from V b into two 
pieces: tree-level t — > W + b+2 jets plus dipole subtraction 
terms, and one-loop t — > W + b+l jet plus integrated 

(2) 

dipole terms, . Together with the NNLO corrections 

to T Al denoted by T ( f\ they add up to the full NNLO 
QCD corrections to T t . 

Finally, we note that throughout the calculation in this 
Letter, the strong coupling constant is rcnormalized in 
the modified MS scheme |19| . And renormalization of 
masses, wave functions, as well as electroweak coupling 
constant are carried out in the on-shell scheme [20(. It 
should be pointed out that the method used here to calcu- 
late the NNLO corrections is similar to the qr subtraction 
method of Catani and Grazzini [l(|. In fact, they both 
employ the universality of infrared divergences and the 
knowledge of rcsummation to facilitate the calculation. 
Total width. For top quark SM decay, the total decay 
width in Gf parametrization scheme [20( at LO is given 
by 



,(0) 



l-3(^f ) 2 +2(^f ) 3 



assuming CKM matrix element | Vtb\ = 1 and mi, = 0. We 
choose m w = 80.385 GeV, G F = 1.16638 x 1CT 5 GeV -2 , 
and TOt 



173.5 GeV 21 



unless specified. Other con- 
stants used in followed calculations include mz, a s (mz), 
and TOb, which are also chosen as in Ref. [Hj]. Correc- 
tions to the LO width considered here include finite b 
quark mass and W boson width effects, S b j and <5j^, NLO 
electroweak corrections, 5ew, NLO and NNLO QCD cor- 
rections, 5q£ D and Sqq D , which are defined as 



r| 0) (l + (5 



>EW 



+ r) {2) 



J QCD t "QCD/' 

where T t is the corrected total width. In Table U we show 
the LO total width together with all the corrections in 
percentage (%) for different top quark mass values. The 
renormalization scale is set to top quark mass. Our re- 
sults agree with the ones shown in previous literatures 
for finite width and mass effects [l| , electroweak correc- 
tions dUlli], NLO QCD corrections Q with the updated 
input parameters. Especially, although using quite differ- 
ent method, our NNLO QCD corrections agree with the 
results in Ref. [f| within the range of the uncertainties of 
numerical calculation, which are of the order 10 -4 . All 
the corrections are stable with respect to the top quark 
mass. 

As mentioned earlier, the NNLO QCD corrections can 

(2) 

be divided into three pieces, T\ with i = 1,2,3. Each 
of them depends strongly on the cutoff parameter tq up 
to the fourth power of In To. While their sum should 
only have weak dependencies proportional to To, and ap- 
proach the genuine NNLO QCD corrections when tq is 
small enough. Thus in Fig. [T] we show the separate con- 
tributions to the NNLO corrections. When tq varies from 
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m t 
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d f d f d E w uqcd °qcd 


172.5 


1.4806 


-0.26 -1.49 1.68 -8.58 -2.09 


173.5 


1.5109 


-0.26 -1.49 1.69 -8.58 -2.09 


174.5 


1.5415 


-0.25 -1.48 1.69 -8.58 -2.09 



TABLE I. Top quark total width at LO and corrections in 
percentage (%) from finite W boson width, finite b quark 
mass, and high orders, including NLO in EW couplings, NLO 
and NNLO in QCD coupling. Mass and width are shown in 
unit of GeV. 



1CP 3 to about 10 -6 , the separate contributions can reach 
as large as twice of the LO width, while the sum keeps 
almost unchanged at the value of about 2.1% of the LO 
width. Stability of such a large cancellation proves the 
validity of our NNLO calculation. On the other hand, 
the NLO QCD corrections have an uncertainty of about 
1.6% of the LO width due to the arbitrary choice of renor- 
malization scale as shown in Fig. [5J which comes directly 
from running of QCD coupling constant a s . After adding 
the NNLO QCD corrections, the scale dependence is re- 
duced to about 0.8%, which makes the predictions more 
reliable. 
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FIG. 1. Separate contributions of the NNLO QCD corrections 
and their sum as functions of the cutoff to, normalized to the 
LO width. 
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FIG. 2. Renormalization scale dependence of the NLO and 
NLO+NNLO QCD corrections, normalized to the LO width 
at central scale fj, = m t . 



Differential distributions. Within our framework we 
can calculate the fully differential decay width of top- 
quark semileptonic decay t — >• W + (l + v)b up to NNLO 
in QCD, which is not possible for the method based on 
calculations of top quark self-energy. Precise predictions 
for differential distributions of top quark decay products 
are of great importance, especially for the measurement 
of top quark mass @ and test of the V — A structure of 
tWb charged current Q • Below we will show several final 
state distributions for t — > W + (l + v)b, including all the 
corrections as in the total width results. We use e + e~ fc^ 
algorithm [23| at parton level with jet resolution thresh- 
old y cu t =0.1 for jet clustering, which is more suitable 
for presentation of the results in top quark rest frame 
as compared to the jet algorithms used at the LHC. The 
shape measurements are more relevant for the experimen- 
tal studies, having both small experimental and theoret- 
ical uncertaities. Thus all the distrbutions shown below 
are normalized to unit area for comparision. For each dis- 
tribution we show results for several cases, i.e., pure LO 
prediction (denoted by LOl), LO predictions plus correc- 
tions from finite mj,, W boson width and NLO EW effects 
(L02), L02 with NLO QCD corrections in addition, and 
further including NNLO QCD corrections. We further 
checked that the NNLO corrections to the distributions 
are also stable against the cutoff tq . 
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FIG. 3. Energy distribution of the charged lepton from top 
quark decay in top quark rest frame. 

From Fig. [3] to El we present the charged lepton 
energy distribution, invariant mass distribution of the 
charged lepton and the hardest jet in energy, in top quark 
rest frame, and two angular distribution of cos(#*) and 
cos(9ij). All of them are normalized to unit area. 9* 
are defined in W boson rest frame as the angle between 
charged lepton and the opposite of top quark direction, 
and 9ij is the angle between charged lepton and hard- 
est jet in top quark rest frame. In each figure the upper 
panel shows the normalized distribution while the lower 
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panel gives their ratios with respect to the one of LOl. 
As we can see, the differences between LOl and L02 are 
small in general, especially for the central region of each 
plot. Both the NLO and NNLO QCD corrections push 
the energy and invariant mass distributions into the cen- 
tral region since the recoil constituents are then massive. 
The NNLO corrections here are about one-fourth of the 
NLO ones, similar to the results of total width. Inclu- 
sive angular distribution of cos(#*) reflects the W boson 
helicity fractions in top quark decay, which can be also 
predicted up to NNLO in QCD through top quark self- 
energy calculations [24j. cos((9*) distributon has been ex- 
tensively studied at both the Tevatron and LHC for test- 
ing potential anomalous tWb couplings induced by new 
physics ■ By a least x 2 fit we get the W boson helicity 
fractions ratio as T L : F + : T- = 0.689 : 0.0017 : 0.309 
using the cos(#*) distribution. The results incorporate 
finite b quark mass and W boson width effects, one-loop 
EW corrections, and QCD corrections up to NNLO, and 
are in very good agreement with the one shown in [24| . 
Our calculations are more helpful for the corresponding 
measurements since experimentalists can include precise 
corrections for the acceptance in different kinematic re- 
gions using our results. As for cos(#y) distribution, QCD 
corrections are more pronounced there since changes of 
the energy spectrum also modify the distribution. 
Conclusions. We have presented the NNLO QCD cor- 
rections to top quark total decay width, which do not 
depend on expansion in W boson mass, and fully differ- 
ential distributions of t -> W + {l + v)b based on SCET. 
One-loop EW corrections as well as effects from finite b 
quark mass and W boson width are also included. All 
together they constitute the current most precise pre- 
dictions for top quark decay, which are helpful for top 
quark mass measurement and testing of weak charged 
current structure. We have implemented the calculation 
into an efficient parton level Monte Carlo program [25j |. 
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FIG. 6. Angular distribution of the charged lepton from top 
quark decay in top quark rest frame. 



in which arbitrary infrared-safe cut can be imposed on 
the final state. Our calculations are complementary to 
the NNLO QCD predictions for top quark pair produc- 
tion Moreover, our method can be widely used in 
studies of heavy-to-light quark decay, including B meson 
semileptonic decay, which will be presented elsewhere. 
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